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Abstract. We consider the tomography problem of recovering a covector field on a simple Rie- 
mannian manifold based on its weighted Doppler transformation over a family of curves V. This 
is a generalization of the attenuated Doppler transform. Uniqueness is proven for a generic set of 
weights and families of curves under a condition on the weight function. This condition is satisfied 
in particular if the weight function is never zero, and its derivatives along the curves in F is never 
zero. 

1. Introduction 

The Doppler transform of a (compactly supported) covector field / = {/j} in R" is given by 

X 9 (x, j fj{x + te)9^ dt. (1) 

Using the Fourier transform, it can be easily seen that one can recover only / up to a differential 
of a function d(j) so that (j) = for large x, see e.g., [13]. In particular, such d(j) (called potential 
covector fields, or simply potential fields) always belong to the kernel of the transform, and this 
follows simply by the fundamental theorem of calculus. The fields orthogonal to potential fields are 
called solenoidal and are characterized by the fact that they are divergence free. In fact, unique 
recovery of the solenoidal part of / is possible in dimension n = 3 from the Doppler transform along 
lines parallel to only two separate planes, although stable recovery in 3 dimensions requires the 
Doppler transform along lines parallel to three different planes, see [12]. Other partial data results 
and methods for recovery of the solenoidal part of / from its Doppler transform can be found in 
[U, [5j, tlOj, and [ilj. 

The attenuated Doppler transform is defined in a similar way: 

i^f{x, 6)= [ w{x + 19, e)fj{x + te)e^ dt, (2) 
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where the weight w is the total attenuation along the ray [x,x — ooO] with an attenuation factor 
aix) 

= e-^--"(^+*^)'^^ (3) 

It came as a surprise to find out in p] (see also [8]) that if a is smooth enough and cr > on supp/, 
then one can recover / uniquely from 1^^,; and there are even explicit recovery formulas of the type 
known for the scalar attenuated transform [9j. 

In this work, we study a more general version of this problem. First, we work on an n-dimensional 
manifold with boundary M diffeomorphic to a ball (where the support of / lies), and instead of 
straight lines, we consider integrals over a general family of curves which we will call T. We can 
choose and fix a global coordinate system on M. The set F will have the following properties. For 
any (x, ^) E TM\{0} there exists a unique directed curve in F through x in the direction of ^, there 
is a unique curve in F through any given pair of points, smoothly depending on the endpoints, and 
dM is strictly convex w.r.t. F. In particular, the curves in F solve the ODE 7 = ^(7,7) with some 
generator G{x,(,). An example are the geodesies of a simple Riemannian metric g on M, see the 
definition below. Then we study the weighted Doppler transform 

/»/(7) = I w{-f{t),j{t))f,{^m^{t)dt, (4) 

and ask the following question: What condition on the weight w{x, ^) and on the family F would 
guarantee that /^u is injective? Clearly, if w = const., then any d4> as above is in the kernel. 
Moreover, if w = const, along any curve in F, but maybe not a constant overall, this is still true. 
On the other hand, that is not true for general w's and there are no obvious elements of the 
nullspace. One consequence of our main result is that if w is never zero, and the derivative of w 
along each 7 G F is non-zero, then the kernel of ly^ is finite dimensional; generically, the kernel 
is trivial and there is a natural stability estimate. In particular, this applies to the attenuated 
Doppler transform on a simple Riemannian manifold, see section HI Actually our result is more 
general than this. We use methods developed in [TTl [6]. 



2. Main results 

Let M be as above; since we have assumed that M is diffeomorphic to a ball we can always work 
in a set of global coordinates on M, which we will do. Let G*(x,i^) for i = 1, ... ,n be a collection 
of smooth real-valued functions on TM \ {0}. We loosely think of G(x,^) = (G*(x,^)) as a vector 
field, although it does not satisfy the transformation law required to actually make it a vector 
field. We will define F as the curves solving 7 = ^(7,7) with given initial conditions on dM. It is 
very convenient to extend M a bit to a manifold with boundary Mi with M™* 3) M, so that the 
curves extend smoothly to Afi (which can always be done if G is smooth up to dM). We can also 
assume that the so extended curves in Mi still satisfy the assumptions in the introduction since 
those conditions are stable under small perturbations, and that Mi is diffeomorphic to a slightly 
larger Euclidean ball. So we now have smooth functions G^{x,^) on TMi \ {0}. We will take this 
as a basis of our definition and continue to describe how the curves in F may be parametrized. 

Fix any smooth Riemannian metric g on M and extend it to Mi. We can always take g to be 
Euclidean, and that is the natural choice if the curves in F are not related to some metric. If there 
is (another) metric on M, and F consists of geodesies of that metric, it will be more convenient to 
make g equal to this metric, but it is not necessary. The raising and lowering of indices, and the 
norms below, are with respect to the metric g. Also, we will often refer to the unit sphere bundle 
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SMi = {{x,e) G TMi I \9\g = 1}, and define 

d±SMi = {(x, 9) G dSMi- ±{v, d)g > 0} , (5) 

where v is tlie exterior unit normal to dMi, and dSMi = {(x, 9) € TMi; x € 5Mi, |^| = 1}. 
Fix a positive function A on d^SMi. Given (x, 9) € d^SMi, we define 7^^£) as the solution of 

7 = G(7,7), 7(0) =x, 7(O) = A(x,e)0, (6) 

defined over the largest possible interval [0,t{x,9)], t{x,9) > 0. Assume that t{x,9) is always 
finite, and also note that there is an implicit assumption in the previous sentence that the maximal 
interval of definition for each of the curves is closed. We now consider T to be the set of curves 
defined by ^ as {x,9) ranges over d-SMi, and define a smooth structure on T by declaring the 
map (x, 9) e to be a diffeomorphism. The generator G defines a vector field G on TMi \ {0} 

by 

Then F consists of projections to Mi of integral curves for this vector field. 

We assume that 7 is uniformly bounded away from zero for every 7 E F as defined in the 
previous paragraph. Then the flow {(7,7)} of G with initial conditions as in ^ deflnes a smooth 
codimension one surface G in TM™*. Assume that Q can be parametrized by (7, 7/|7|g) G 5* Mi, 
i.e., V(y,a;) € 5'Mi, there are unique t>0, {x,9) € d-SMi, smoothly depending on {y,uj), so that 
y — 7(*)) ^ — ii't)/\i{'t)\g with 7 = ^x,e- Then X{y,uj) = \'^\g is a well defined smooth positive 
function on S'Mi which can also be extended to TMi \ {0} by homogeneity of order —1. That is, 
X{y,9) = X{y,9 /\9\g)/\9\g. This is an appropriate extension because with it X{x,9)9 is the velocity 
vector for some curve in F for every (x, 9) € TM. We also extend the definition of ^^fi as in ([6]) 
but now X does not need to be on dM. 

Given any (x, 9) G TMi \ {0} let the domain of the maximally extended curve ^xfi be given by 
[r_(x, r+(x, 0)]. We say that (G, A) deflnes a simple system, if for any x € Mi the maps 

{{t,9)- 0<t< T+{x,9), \9\ = l}3 {t,9) ^ expx,+ {t,9) = 7^,e(t) G Mi \ {x} 

and 

{{t, 9); 0<-t< -r_(x, 9), \9\ = 1} 3 {t, 9) ^ exp,__(t, 9) = jx^g{t) G Mi \ {x} 

are diffeomorphisms, depending smoothly on x. For each x, we will denote the maps given above, 
taken together on the domain including t = 0, by exp^{t, 9). 

An example of G is the generator of the geodesic flow corresponding to g on Mi (or to any other 
metric g on Mi). Then we can choose A = 1 (or A = |^|g/|^|§). The simplicity condition is fulfilled 
if g is a simple metric \1?>\ I17j. Another class of examples are magnetic systems on {Mi,g), where 
the geodesic equation is perturbed by a term representing a Lorentzian force, always normal to the 
velocity [3]. Then the fiow preserves the energy level, so we can choose A = 1 as well. 

As stated in the introduction, the main object of our study will be the weighted Doppler transform 

IgXw ■■ L\M) ^ L\d.SA'h,det{g) {9,u)gdSxdSg) 

given by the formula Here and in what follows L?'{M) refers to the space of covector fields. 
Also, the family of curves F has been identified with d-SMi, and the measure det{g) {9, v)gdSxdSe 
is chosen to provide a convenient structure on F (other measures could also work). The measure 
det((7) {9, v)gdSxdQ is invariantly defined but depends on the choice of g. The fact that the weighted 
Doppler transform has this mapping property can be established by considering the operator N = 
Iq ^ y,lG,\,w which is a ^J}0 of order —1. 

Before stating our main results we require one more definition. 
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Definition 1. We say that w satisfies the elliptic condition, if there exists an open set U C Q 
satisfying 

{C e T*M \3eeU with {C, 9) = 0} = T*M, w\u + 0, (8) 
and for any (x,^) G T*M, there is no covector h £ T*M such that 

Glogw{x, 9) = hje^ (9) 

for all 9 normal to ^ with (x, 9) G U . 

Condition ([5]), combined with the non-simpUcity condition, is the necessary and sufficient condition, 
provided that U is open, for the ehipticity of Iq ^ w^G,x,w when acting on functions, see [6]. 

The case of dimension n = 2 is special because then there are only two vectors 9 such that 
{x,9) £ Q is normal to any given ^. If U were a proper subset of Q, then for some ^ the set of 
(x, 9) £ U with 9 normal to ^ would contain at most one vector, and therefore ([9]) is trivially true 
for some h. Thus we observe that if w satisfies the elliptic condition and the dimension is 2, then 
it must be the case that U = Q. Taking this into consideration, in dimension 2 we see that the 
elliptic condition merely states that for any (x, 9) G TM \ {0} 

Glogwix,Xix,9)9) + -^23^^G,\ogw{x-\{x-9)9). 

A(x, —9) 

In particular, if A(x, 0) is even in 0, which is true in the Riemannian case (where A = 1), then the 
condition above simply states that Glogt(;(x,0) should not be an odd function of the 9 variable. 
Note that the following condition is sufficient to imply that w satisfies the elliptic condition 

u; / 0, Gty 7^ 0. (10) 

This says that w satisfies the elliptic condition if it does not vanish, and its derivatives along the 
curves in F also do not vanish. In the case of the attenuated Doppler transform on a Riemannian 
manifold where w is given by a Riemannian version of ([3]), (jlOp is satisfied because then Qsw = 
—a{x)w{x,9), where o" > is the attenuation. Notice that a > can depend on the direction as 
well, and then (jlOp still holds. For more details, we refer to section HI 

On the other hand, suppose that /i(x) is any smooth covector field on Mi. Then ([9]) can be 
solved for logw{x,9), and therefore for w{x,9), by integrating along curves (7,7) C TMi where 
7 G F. For any ip G C°°(Mi) that is zero on dMi we then have 



lG,x,Ui^h + dij){j) = I u;(7(t),7(i))(V'(7(i))/i(7(t))i +dV(7(*))i)7(iFdf 

wil^ hj 'y{ty — ip Gw dt 
wip hj j{ty — ipw hj j{ty dt = 0. 
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For the second equality we used integration by parts while the third follows from ([9]). Thus we see 
that some condition on w is required for injectivity. In fact, with the elliptic condition we have our 
two main results. 

Theorem 1. Let w satisfy the elliptic condition, and let G, A, and w be real analytic on a neigh- 
borhood of the curves {'^xfiiixfi) for every {x,9) G U where U is the set from DefinitionUi Also 
suppose that (G, A) forms a simple system. Then Ig,x,w is injective. Further, if lG,\,w{f){lxfi) = 
for all (x, 9) £U then / = 0. 
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Theorem 2. Let A G C°°(S'Mi), G G C°°(TMi \ {0}) /orm a simple system. Let w G C°°{TMi \ 
{0}) satisfy the elliptic condition. Then the kernel of Ig,\,w is finite dimensional. Moreover, if 
Ig,x,w is injective then we have the following stability estimate 

WniLHM) + MlHM) < Cm^-^Jf,0]\\uHM,) for any f G L\M) (11) 

and for any (G, A, w) in a neighborhood of (G, A, w) with some constant C > 0. Here and 
<j) are the solenoidal decomposition of f on M defined below, and N is the normal operator also 
defined below. 

The normal operator N appearing on the right hand side of (jlip is equivalent to a continuous 
operator applied to Iq ^^^^ really can be called a stability estimate even though it does 

not provide such for Iq xw directly. Theorem [2] shows in particular that if Ig,x,w is injective for 
some {G,X,w) G C°°, then for every {G,X,w) in some neighborhood of {G,X,w), Iqxw ^^^^^ 
injective. Combined with Theorem [H this gives uniqueness for an open and dense set of {G,\,w) 
with respect to the G^ norm, which we call generic uniqueness. It is also true that the C°°'s in 
Theorem [2] may be replaced by G^ for A: ^ 4. 

3. Proofs 

We start with the Hodge type decomposition of vector fields in M, already used in the statement 
of Theorem [2j Given / G L^(M), we write 

f = r + dct>, (12) 

where 5f^ = and (f) G Hq(M). Here 5f = V^fi, where V is the covariant derivative. The function 
(j) solves the elliptic problem 

AcP = 5f, (PIqm = 0, (13) 

where A is the Laplace-Beltrami operator associated with g. Then G L'^(M), and (p G Hq{M) 
by elliptic regularity, see pO]. We will also need to use the solenoidal decomposition of / on Mi, 
which can be obtained by extending / as zero on Mi\M and replacing dM by dMi in (I13p . In fact, 
we will generally consider any function space on M as the subspace of the corresponding function 
space on Mi where the functions are are extended as zero on Mi \ M. 
Plugging (fT2]) into the definition (jl]) of Ig,\,w we get 



Ig,a,-/(7) = j wf'^^Ut + j w^j^dt 
wf^j^dt- j (j)Gwdt. 



(14) 



The problem therefore can be reformulated as follows. Given a covector field / and a function (p, 
consider the pair [/,(/>]. We say that [f,4>] is a solenoidal pair if 5f = 0. We then study the ray 
transform 



![/,</'] (7) = J {wfjj^ +a^)dt, (15) 

where w{x, ^) and a{x, £,) are given functions on TM. This is actually an integral of w{x, £,)fj{^)C'' + 
a{x,^)<j){x) along the maximal integral curves of G. We are mostly interested in the special case 
when a = —Gw, and in this case the second equality in (fH|) shows that I[/ + d4>, ip] = I[f, + Tp]. 
However we can also consider the more general case when a may not be related to the weight w. 
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Remark 1. If we reformulate the elliptic condition in terms of a, then the proofs of propositions [T] 
and [2] will work in the general case when a does not correspond to —Gw. The appropriate refor- 
mulation is to simply replace Glogw(x,^) in Q by a{x,6)/'w{x,6). 

We consider [/, (p] as an element of the space (Mi ) with norm 



where dvg is the volume form associated to g, and |/| is the norm of the covector / in the metric 
g . Solenoidal pairs then are exactly those orthogonal to potential ones given by [dV',0], where 
V' € H^iMi). We will denote the space of solenoidal pairs on Mi (resp. M), which is a closed 
subspace of L^(Mi) (resp. L^(M)), as Lg(Mi) (resp. Lg(M)). We will also need to consider 
Sobolev spaces of pairs [f,(f>], and will denote these, as in the statement of Theorem [21 with bold 
letters. We view I as an operator I : L^(M) — > L^(9_S'Mi, det{g) {9, i')gdSxdSg) , and we set 



As mentioned above, the fact that N is a ^DO of order —1, which is established in Proposition [H 
proves that I does indeed map L^(M) to L'^{d-SMi,det{g) {9,iy)gdSxdSg). When it is necessary 
to distinguish between the operators I and N corresponding to different triples {G, A, w), or just to 
emphasize the dependence, we will sometimes write Ig,x,w and ^g,x,w 

We say that N is elliptic on solenoidal pairs, if diag(A5, 0) and N form an elliptic system of 
operators, where A is any parametrix of A in Mi. In other words, we want A5f = and N[/, </>] = 
to be an elliptic problem. 

Proposition 1. // (G, A) is a simple system, then N is a '^DO of order —1 in the interior of Mi. 
If a = —Gw and w satisfies the elliptic condition, or a and w satisfy the modified elliptic condition 
from remarkUl then N is elliptic on solenoidal pairs. 

Proof. We follow the arguments in [6], where weighted integrals of functions of x are considered. 
Note that we can think of ^-^ in w{x,^)fj{x)^^ = {'w{x,(,)^^) fj as part of the weight. Then we can 
apply the analysis in [6] (see the proof of Proposition 2 there) to w^^'^ fj^ for each jo, and also to 
a(/>. Summing the results we thus obtain the symbols of the corresponding \I'DOs. On the other 
hand, N is an operator acting on pairs, unlike [6]. The latter case is done in [3], where F consists of 
magnetic geodesies; see Propositions 3.9 and 4.1 there. Thus we get that N is an integral operator 
of the form 




N = ri. 




(16) 



with 




(17) 
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Here is an elliptic factor which is given by the Jacobian of the change of the variables {z, uj, t) i— > 
{x,6), where {z,uj) € d-SMi, ^z,uj{t) = x, 'yz,Lu{t) = A6'. If F consists of geodesies, or even magnetic 
geodesies as in [3], then = 1 by the Liouville theorem. 

Following OE], we deduce that N is a ^'DO of order —1, and for the principal symbols ap{Nij) 
we get 





1 


\\x.. 


, 9)9''9'\w{x, X{x, e)e)\^J^{x, e)d{{C, e)) da^iO), 


ap{Niof{x,0 = 


j 


X{x,t 


9)9'' w{x, A(x, 9)9)a{x, A(x, 9)6) j\x, 9)6{{^, 9)) da^{9) 


ap{Noiy{x,0 = 


/ 


X{x,( 


9)9'a{x, A(x, 9)9)w{x, A(x, 9)9)j\x, 9)6{{^, 9)) da^{9), 


crp{Noo)ix,0 = 


/ 

J So: Ml 


\a{x, 


X{x,9)9)\^J^{x,9)6{{C,9))da^{9). 



It follows that 



f,X{9)9^+ '^y:)Zn. <P = 0- (20) 



K(N)[/, </.],[/,(/>]) = / \X{x,9)w{x,X{x,9)9)fj{x)9^+a{x,X{x,9)9)cl>{x)\'j\x,9)6{{t9))da,,{9). 

Js^Mi 

(18) 

Note that in the case of dimension n = 2, (jlSp is actually just a sum of two terms. The principal 
symbol of A6 is |'^|~^'^*. Fix x and 7^ 0. Assume now that 

he = (19) 

and (lisp vanishes. By ([8]), the set of unit vectors 9 normal to ^ so that (x, A(x, 9)9) G U, is an open 
subset of {\9\ = 1, ^) = 0}. In particular, w{x,9) ^ for all such 9. So we have 

a{X(9)9) 
w{X{9)9)' 

Since x is fixed, we suppressed the dependence on x. If cp = 0, then (fTUj) and (pUj) easily imply 
/ = as weh, and we are done. If 7^ 0, then we get a{X{9)9)/w{X{9)9) = -cj)-^ fjX{9)9^ for 9 
as above. Depending on which case we consider either this last equation contradicts the modified 
elliptic condition given in remark[Tl or a = —Gw and then since Glogt(;(x,0) = {Gw)/w = —a/w 
we have a contradiction of the original elliptic condition. □ 

Proposition 2. Under the assumptions of Proposition{l\ the following a priori estimate holds for 
any solenoidal pair [/, i?5>] G L^(M). 

II/IIl2(m) + II0IIl2(m) < C (l|N[/, (/']||hi(Mi) 
Proof. Without loss of generality, we can assume first that /, (j) are smooth. Also, let Mx/2 and 
M3/4 be other manifolds with smooth boundary so that M (s M*"* <e ^t^A ^ -^i""*- 

Since the pair of operators (N, diag(A(5, 0)) is elliptic, there exists a parametrix for the system 
on the interior of Mi. We will refer to this parametrix, which is a system of ^'DO's, as {A,B). 
Thus we have for any pair \h,ijj\ € L^(Mi) 

^N[/i,'0] +5[A(5/i,0] = [h,^lj]+K[h,i}] 

where K \s a smoothing operator. 

Now, take a cut-off function x ^ Co°(M{"*) that is equal to 1 on M3/4, and let /|^^ be the 
solenoidal projection of / (extended as zero to Mi \M) in Mi . Then since supp((5x/|fj ) C M\M3/4, 
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by the pseudolocal property of ^fDO's, B[A6xfM-^^, 0] smooth on a neighborhood of Mi/2- Thus, 
applying (A, B) to (N, diag(AJ, 0))[/|.j^ , (/>] we have 

xAnxfk^<l>] = [/m„0] +i^[/M,,'/'] in Ml/2, (21) 

where K has a smooth properly supported (in Mi) kernel. 
Recall that = / — dcpMi, where 

A^M^=Sf in Ml, cj)M, e H^{Mi). (22) 

Thus the term K[f^,j^,(j)] from (j2ip can be replaced by K[f,(j)] with a different -ftT with the same 
property. 

We will write (A^^)^^ for the solution operator of the Diricihlet realization of the Laplace- 
Beltrami operator A on Mi so that with this notation {A'^j_^)^^6f = (pMi- Now, let A~^ be a 
parametrix for the Laplacian on Mi so that {{A^j^)^^ — Ap^)6 f = Kif on M1/2, where Ki is another 
operator with a smooth properly supported kernel. Since /|,^^ = / — d</>Mi = f — dKif + dA~^6f 
and xd = dx — (dx) we get that 

AN[xfM,A] = ^N[/, </>] +^N[dxKi/, 0] -^N[(dx)ifi/, 0] +AN[dxAp'6f, 0]-AN[{dx)A^'6f, 0]. 

By (lisp , the principal symbol of N vanishes on potential pairs, i.e., pairs of the form [d'(/',0], and 
so all of the terms on the right hand side except the first may be grouped together and written as a 
single system ^fDOs of order —1 applied to [/, </>]. Applying all of these comments to ()2ip we have 

AN[f,<p] = [flj^,<P]+K[fA] in Ml/2 (23) 

where is a new ^DO of order —1. 

Now, denote the projection [/, i?i>] 1— > / by vri. From (j23p we have 

dcf>M, = vri ( - ^N[/, </.] + [/, 0] + K[f, <P]) in M1/2. (24) 

In particular, 

dcPM,=7Ti{-AN[f,<P]+K[f,4>]) in Ml/2 \M (25) 
since [/, </>] = outside of M. We integrate the equality above to get 

(PmAx)= f\M^^if^^]+K[f,<P]),c,)ds + K2f, xeMy2\M, (26) 
Jo 

where [0,1] 3 s Cx{s) is any curve outside M connecting x and a point on dMi/2, smoothly 
depending on x. If dMi/2 is at a fixed distance to dM, which we can always assume, then Cx can 
be chosen to be the normal geodesic to dM. The operator K2, given by K2f{x) = 0Mi(cx(O)), 
is smoothing since (j)Mi is harmonic on a neighborhood of 5Mi/2. In particular, K2 ■ L?'{M) — > 
H\Mi/2\M). 

Equalities dM]), dM]) yield 

H-i(Mi) ) • (27) 

Using (f27|) and the trace theorem we can also estimate \\4>Mx \\H^/'^(dM)i ^nd since / is solenoidal on 
M, (j)Mi solves A(/)Mi = in M. Using the standard estimate for solutions to the Dirichlet problem 
we thus have 

UMAn^i^M) < C (||N[/,<A]||hi(a./o) + \\[f^'P]\\n-HM,)) ■ 
Combine this, ([23|) and to get the estimate in the proposition. 

□ 
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Remark 2. The arguments of the previous two propositions will still work for (G, A, w) € if k is 
chosen large enough. When working with the ^'DO calculus we actually only need to expand the 
symbols up to a certain order assuming that the remainder is still sufficiently regular. 

Proof of Theorem\M First take any solenoidal pair [/, (f)] S Ker(NG'^;^_^) nhKM). By Proposition [2] 
we have that 

\\[f,mL^(M)<C\\[f,mu-HM,y 

Since the inclusion i : L^(M) —>■ Ii^^{Mi) is compact, this implies that the intersection Ker(NG^A,i«)n 
L^(M) is finite dimensional. Using (fn|) and the fact that Ker{lG,x,w) = Ker(NG_A,w)) we see that 
KeT{lG,x,w) is also finite dimensional, which establishes the first statement of Theorem [2j We also 
make the observation that Ker(NG,A,t«)nL2(M) C C~(Mi) (C°°(Mi) is simply the space of smooth 
pairs [/,(/>]). We can see this by using a parametrix for the pair of operators (N^^^, diag(A(5, 0)) 
For the stability estimate we will make use of the following lemma which is similar to [19^ 
Prop. V.3.1]. See also [161 Lemma 2]. 

Lemma 1. If X , Y, and Z are all Banach spaces, A : X ^ Y is a continuous and injective linear 
operator, K : X ^ Z is a compact linear operator, and we have the estimate 

\\x\\x <C{\\Ax\\y + \\Kx\\z) VxeX, (28) 

then in fact we have 

\\x\\x<C\\Ax\\y VxGX 

Proof of Lemma[l\ Suppose that the conclusion were not true and so we could find a sequence 
{xn}'^=i C X with = 1 for all n, and such that — >■ as n ^ oo. By compactness of 

K, there exists a subsequence {xnk}^=i ^^^^ {^^n,,}^=i converges as /c — > oo and is therefore 
Cauchy. Thus using (f28l) and the fact ||ylx„j.|[y ^ we see that {xn^.}^]^ is Cauchy. Thus {xn^} 
converges to some x (z X. The continuity and injectivity of A show that x = 0, but this is a 
contradiction of the original assumption that ||x„|[x = 1 for all n. □ 

To apply the lemma in this situation first note that since Ig,\,w is injective and Ker(NG_A,M;) H 
L2(M) C C~(Mi), by dm) if [/^ (/.] G Ker(NG,A,«;)nL2(M) then f'+d^ = 0. Since f is solenoidal 
in M this implies A(p = in M, but also (/'Igjv/ ~ so = which in turn implies = 0. 
Thus N^A.u, is injective on L2(M) and so we can take X = L2(M), Y = Hi(Mi), Z = H-i(Mi), 
A = l^G,x,wi and K the inclusion map from X to Z. Using the estimate from Proposition [21 we 
establish 

\\[r,mL^(M) < cmcxAr^mn^iM,) v[r,<^] g lum) (29) 

for some constant C > 0. If / G Lp'{M), we may now take its solenoidal decomposition f = f^ + d(f) 
on M. Applying ^ to [f'A] yields 

wrwi^iM) + uwlhm) < c\\NG,xAr,mHHM,y m 

If we note that by ([Hj), Ng'^;^^^[/, 0] = ^G,x,w[f^ , ^P] we can see that this is the stability estimate 
pT]) in the case (G, A, w) = {G, A, w). 

To complete the proof we argue that for small enough perturbations (G, A, w) of (G, A, w) in G^ 
the estimate ([30]) still holds with a constant G independent of (G, A, w). Note that in the case that 
(G, A, w) are only G^, we may still define Ig,x,w and the normal operator ;^ which is still given 
by (jl6p . We will follow [6] and establish the following proposition. 
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Proposition 3. // {G, A, w) G C°° is as in the statement of Theorem 0, and (G, A, w) £ 

with \\{G,X^w) — {G,X,w)\\qs < 5, then for all 6 sufficiently small (G, A) is still a simple system 
(where the diffeomorphisms are only C^j, and there exists a constant C > 0, which depends only 
on {G, A, w) and the dimension n such that 

||(Ng,a,--N^^3,_J[/,<A]||hi(aa) <C'<^(II/IIl2(a./) + II</'IIl2(m)) V/,<AgL'(M). (31) 

Before proving Proposition [3] we note that the proof of Theorem [2] follows if we take a new 
constant C which is the product of the two constants from (I30p and ()3ip and then make 6 < 1/2C. 
Using both (f30l) and (f3T]l then gives 



1 

2 



< c \\Ng~,jr,mHHM,) + hn L^{M) + \mL^{M)) 



and so 

WfhHM) + Uh^m < 2C ||N^,x,^[/,0]||hi(Mi) V/ G L\M) 

and for all (G, A, w) in the 6 neighborhood of (G, A, w) with respect to the G^ norm. This completes 
the proof of Theorem [2j □ 

Proof of Proposition O The proof will consist of a careful comparison of the integral kernels corre- 
sponding to ^G,\,w and ^ ^ given by (fT7|l . We begin by considering the two maps 

: ^ (r,u;) = fsign(t) |exp,(t,e) - x|,sign(t) ^^H^Lgl^l^) . (32) 



Here exp^ may correspond to either (G,A) or (G,A), and the corresponding maps given by ([32]) 
are denoted by and respectively. We extend both (G, A) and (G, A) sufficiently far, possibly 
beyond Mi, so that these two maps have the same domain. The {t,6) coordinates are loosely 
speaking polar coordinates with respect to the family of curves generated by G (or G) while (r, lo) 
are actual (Euclidean) polar coordinates centered at x. 

Our first order of business will be to estimate F~^ — F~^ in the Gj ^ ^ norm on the set of (x, r, uj) 
such that X + rco G Mi. We will accomplish this in three steps. First we will estimate {F^ — F^), 
and then using this we estimate {F~^ o F^ — Id). Finally, we precompose by to obtain the 
desired estimate. 

For the first step we use the following standard estimate from ODE theory which is proven in 
[2] and also applied in the same manner in [6J. 

Lemma 2. Let x and x solve the ODE systems 

x' = F{t,x), x=F{t,x), 

where F, F are continuous functions from [0, T] x U to a Banach space B, where U d B is open. 
Let F he Lipschitz w.r.t. x with a Lipschitz constant k > 0. Assume that 

\\F{t,x) - F{t,x)\\ <6, Vt e [0,r], Vx G U, 

and that x{t), x{t) stay in U for < t < T . Then for <t <T 



\xit)-xit)\\ <e^'*||x(0)-x(0)|| + ^( 



e^'-l 
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We apply the Lemma [2] to the systems 

7'(i) = m fit) = m 

considered with the same initial conditions 7(0) = 7(0) = x and ^(0) = ^(0) = \{x,6)6 to obtain 

||7x,e - IxfiWcl , , + Wixfi - IxfiWcl , , = 0{5). (34) 

To get the estimates for the x and derivatives we first differentiate the systems ()33p with respect 
to the initial conditions, and then apply Lemma[2]as well as the hypotheses about (G, A) and (G, A). 
For the t derivatives we appeal directly to the original system (j33|) and use the assumptions on G 
and G. For the mixed derivatives we use a combination of these techniques. 

At this point note that (i34|) proves the statement that (G, A) is still a simple system, and together 
with the fact that we may rewrite (j32p in terms of ^^fi (resp. ^x,b) it also shows that 

\\Fx{t,e) - Fx{t,e)\\ci =0{5). (35) 

Now note that {F^^^ o F^ — Id) = {F~^ o Fx — F~^ o F^), and so, working in some appropriate set 
of local coordinates for uj S S""-*^, we have 

(F-i o Fx - Id)it, 9) = (^^' DF-\sFx{t, e) + {l- s)Fx{t, 9)) ds^ ■ {Fx{t, 9) - Fx{t, 9)). 

Taking derivatives of his last equation we see that \\F~^ o Fx — Id\\(j3 ^ can be bounded in terms 
of §5i), and \\DF-^\\c3 . By ([35]), ||^;|[c;3 and \\DF-^\\ (j2 are uniformly bounded for (G, A) 

x,r.iij x,t,9 x,r,uj 

in a sufficiently small G^ neighborhood of (G, A). Thus, if we precompose F~^ o Fx — Id with F~^ 
we see that 

|IF-I(r,^) -F-i(r,.;)|b3_ = 0{S). (36) 

Now we use the maps F~^ and F~^ to change variables in ([T71) for N = 1^g,x,w and N = ^ ^ 
respectively. We will call the Jacobian determinants of these changes, which depend on x, J{r, uj, x) 
and J{r,uj,x). By ([36]) we have that \\J{x,r,u;) — J{x,r,u;)\\(j2 = 0{S). We also introduce a 
cutoff function x ^ C^{Mi) equal to 1 on M, and since the integrands in (jl7p are zero outside of 
M we may multiply them by x(3^ + ''<-^) without changing the values of the integrals. Thus the first 
equation for N from ()17p becomes 



{Ni^f)'ix)= xix + ru;)A\\{x,r,u;)Mx + ru;)drdu;. (37) 

where 

^11 (^) ~ 

( A(x, 9)w{x, A(x, 9)9)w{'^xe{t),ix e{t))0'% e{ty' j\x, 9)J{x, r, uj)\ 

We have the same equation for N when tildes are added in appropriate places. Also, we can find 
similar equations for A'^io, Nqi, and A'^oo, and the following analysis proceeds in the same manner 
for those operators. We will only examine A^n in detail. 

By changing variables (r, a;) 1-^ (— r, —uj) we see that if Af^ is odd with respect to (r, uj), then ([37|) 
is always zero. Thus we may replace x^u by (^e)ii = xi^ + ''"^){^ii{x,r,uj) + Ai\(x, — r, —uj))/2 
which is even with respect to (r, uj) . Having done this we take the linear approximation in r for 

{Ae)f[ 

{Ae)f[{x,r,u;) = {A,)f[{x,0,uj) + rR^'' {x,r,uj). (38) 
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Provided that A^{[ &Cl^^ we see that i?**' G C].^^. Changing to Cartesian coordinates y = x + ruj 
in dSZD yields 

JR" \y -x\ \y - ^ 

In the same manner as above we find kernels A^, and R associated to {G,\,w). Thus, comparing 
the two operators we have 

((iVii - Nii)fr\x) = 2^^ ((^e)l^(x,0, ^) - (ie)l^(x,0, ^^^dy 

+ 2 / f \y-xl - i?'^'(x, |y - x|, ^^)') ^^^dy. 

Jr" V \y -x\ \y -x\ J \y- ^ 

(39) 

Let A{x,uj) and R{x,r,u;) denote respectively the portions of the first and second integral kernels 
from pop contained in parentheses and assume for the moment that 

\\Aix,u;)\\ci^^ + \\R{x,r,u;)\\ci^^^^ < C6 (40) 

for a constant C. Assuming / € L?'{M), we may estimate easily the second integral in (I39p and its 
first derivatives in L^(Mi) hy C5 since the kernels of both will have only integrable singularities. 
Similarly, the first integral in (jSOh may be estimated in L^(Mi) by C5. However, the derivatives of 
the first integral pose a problem since when the kernel is differentiated its singularity is no longer 
integrable. However we may still use the Calderon-Zygmund theorem to estimate the first integral 
with a differentiated kernel in L^, and by [3 Theorem XI. 11.1] this is the derivative of the integral. 
Therefore all that remains to prove the proposition is to prove the estimate (j40p . 

To prove (|10]) we first note, from the definition of A[x,uj) and R{x,r,uj), that it is sufficient to 
prove 

\\Af[{x,r,u;)-Af[{x,r,u,)\\c2^^^^=0{6). 

We require since R{x,r,uj) already involves one derivative. It is thus sufficient to show that the 
differences of each of the corresponding terms in the definitions of Af[ {x, r, to) and A^^^ (x, r, cj) (see 
below (j37|) ) are 0{6) in ru- Indeed, this has already been shown for {J{x, r, to) — J{x, r, uj)). The 
other terms all involve {t{x,r,uj),9{x,r,uj)), and so we can estimate them in C^,,,cj by using (f36]l 
together with an estimate of the difference in C^jg, and a C^'^ bound for {G,X,w) (the bound 
must be in C^'^ when we are looking at a = —Gw). Here and in what follows the C^'^ norm is 
the usual ^ q norm plus the Lipschitz norm of all kth derivatives. 

As an example of this method we estimate X{x, 9(x, r, u)) — A(x, 0{x, r, uj)) in ^ w where 9 and 
9 represent the projections onto the 9 component of and respectively. First we have 

\\X{x,9{x,r,uj)) - X{x,9{x,r,uj))\\ci^^^ < \\X{x,e{x,r,uj)) - X{x,9{x,r,uj))\\ci^^^ 

+ \\X{x,e{x,r,uj)) - X{x,9{x,r,uj))\\ci^^^- 

We now see that the first term above can be bounded by (p6]) and a C^'^ bound on A. The second 
term is 0{5) by the hypotheses, and so in fact the entire right hand side is 0{5). The same method, 
also using (j34p . works for the difference between the w and w terms. For the difference between a 
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and d, which is not required for Nn but does arise when we consider A'gi and Nqq, the same method 
works, but because a already involves one derivative of w we need that \\w — w\\(js = 0(6) and 

require a bound on w in Cffg. Finally, for (>-J^) we use ()34p , ()36p , together with the smoothness 
of the boundary of Mi . This completes the proof of Proposition [3l 

□ 

We next move on to the proof of Theorem [TJ Here we will not make use of the normal operator 
N, but rather analyze Ig,\,w itself in a neighborhood of a single curve of F. Our method is analytic 
microlocal in nature and relies on the calculus of analytic wavefront sets which we will denote by 
WFa([/, (j)]) (see [E]). The same technique has previously been applied to other problems: see [6], 
[Is] . The main step is to establish the following proposition. 

Proposition 4. Let G, A, w, and U be as in the hypotheses of TheoremUi Suppose [f^, cp] € L^(Mi) 
is a pair such that is solenoidal in the interior of M, (p]i'yx,d) = for every {x, 9) G U , and 
supp([/'^, C M. Then [^^^(j)] is analytic in the interior of M. 

Proof. This proof will follow closely that of the corresponding results in [6] and [18]. Take any 
(^0) ?'^) £ T*M with xq in the interior of M. We aim to show that 

(rco,e°)^WFA([r,<A]). 

The result then follows since the projection of WFa([/^, (j)]) to M is exactly the set of points where 
is not analytic. 

By the elliptic condition, there exists (xq, ^o) S U normal to (xq, i^)- Since (xq, ^o) is in Q, there 
must be a point {pq^ujq) e d-SMi such that 7po,cjo passes through (xc^o)- We will now define 
coordinates in a neighborhood of 7po^^(j. First, let us take some coordinates (lji, ... ,ujn-i) on Sp^Mi 
centered at ujq. Then by the simplicity assumption (wi, ... ,u;„_i,t) i-^ exppg(t,t(;) is the inverse of 
a coordinate map on a neighborhood of 7po,ajo- We will write x' = (x^, ... jx""-*^) = (w^, ... ,0;"^-'^) 
and x"" = t. Next we translate the coordinates so that xq = and now assume that they are defined 
on the set 

V = {\x'\ < e, /_ < x" < /+} C Ml 

where I- and 1+ are such that (x', 1-) and (x', /+) lie outside of M for |x'| < e. On a small enough 
neighborhood we may assume that these coordinates are analytic because G and A are analytic on 
a neighborhood of ^p^^ujo- Also note that in these coordinates (0, (0, 1)) = (xq, ^o)- 

Additionally, if we fix any 6' = (6^, ... , 9"'~^), then for such 6' sufficiently small the map (z', t) 1— > 
exp(2',o)(^' (^'> 1)) defines another analytic coordinate system, which depends analytically on 0', on 
some subset of V. In particular, if we restrict to \z'\ < 7e/8, then for 6' small enough the curves 
7(2',o),(6»',i) reach points outside of M while still inside V. Thus since U C Q is open, by taking 
\6'\ ^ 1 we can ensure that (/)](7(2; o),(6i',i)) = 0. 

Now we introduce a sequence of cutoff functions xn £ C^{W"'^^) such that supp(xAf) C {\z'\ < 
3e/4}, XNiz') = 1 for \z'\ < e/2, and 

\d^XN{z')\<{CN)\-\ , Vz',|r|<iV. (41) 

It is possible to construct such cutoff functions, see j2H Chapter V, Lemma 1.1]. 

For all ^ in a complex neighborhood of and h > we multiply (/)](7(^/ o),(6i',i)) = by 
XN{z')eh^^ and integrate in z' to obtain 

J I e*^''^^X^(^')(^/;7/,.,o),(eM) + «0)dtd/ = 0. (42) 
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As observed above, {z' , t) provide coordinates on a subset of V, and outside this subset the integrand 
in the following formula is 0. Thus we may change coordinates in (142p to obtain 



eii''(^'nO(^aN{x, 9') // b^ix, 9') + c^(x, 9') cj))dx = (43) 

where for sufficiently small 9' we have that a^v and cjsr are analytic and independent of on a 
neighborhood of 7po,cjo! ^^'^ satisfy (j4ip . with a new constant C, in all variables everywhere. Also, 
aN{0,9') =w{0,X{9',l){9',l)) and Civ(0, 61') = a{0, X{9' ,1){9' ,1)). Finally, the components 6^' (x, 6*') 
of the vector field b are analytic everywhere and we have b{0,9') = A(0, {9', 1)){9', 1). 

Now we will choose 9' based on ^. For convenience we will write 9 = {9', 1). First, without loss 
of generality we may assume that = dx"^^. Now, following [6] and [18] we are able to define an 
analytic mapping 9'{^), and corresponding 9{^), for in a neighborhood of = dx^~^ with the 
following properties. 

9(0-^ = 0, ^"(6 = 1, 0{e) = 9o, 

and 

Note that in order to define a function with all of these properties we need that 9q-(,^ = 0. Replacing 
9' by 9'{S,) in (03]) we arrive at 

J e^'^^''^^\dNix,OffVix,0+CNix,0^)dx = (44) 

where ajy, and cat have the same properties as a^, V , and cat respectively, and 

i^{x,i)=z'{x,9'{i))-i. 

Our next step will be to apply the method of complex stationary phase (see [E]), and so we must 
analyze the critical points of this phase function. In fact, ip is precisely the same as the phase 
considered in [6j, and so we take two results from that paper. Firstly, we have 

^'.^(Ce) =Id, (45) 

and second we have the following lemma. 
Lemma 3. There exists 5 > such that 

^ i^dv^O forx^y, 
for X £ V, \y\ < 6, \^ — < 6, ^ complex. 

It is in the proof of this lemma that 9{(^) • ^ = is used. 

Now for y and 6 as in lemma [3] and r] E B{^^, (5/2), we multiply ()44p by 

where x is the characteristic function of 5(0,(5/2), and then integrate in This yields 

ei'^iV'-'^^O^ANix, v) fj b'ix, + CNix, t v) ^) dx d^ = (46) 



where 



^{x, y, ri, = -V'(y, e) + i^ix, + 2 - 



An and Cn simply incorporate the dependence on rj into cin and cat. Let us first consider the 
portion of the integral in ()46p where — y| > S/C^ for some constant > 1. By lemma [3] there 
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is no real critical point for the function ^ i-^ ^{x, y, r], ^) on this set, and thus we can estimate that 
portion of the integral as follows. 



ei^(y'-'^^^){AN{x, ri) /; P{x, + Cn{x, ^ v) 4>) dx d^ 

\x-y\>5/C^ 



(47) 



{C'^Nhf + C''Ne-'/^^^'^ 

for a new constant C^. To obtain (j47p we integrate by parts N times in ^ using the facts that l*!'^! 
is bounded below on the domain of integration, 

and on the boundary of the domain of integration in (where l^^ — ?7| = <5/2), 

et*(^,s/,^,0 is expo- 
nentially small in 1/h. We also made use of the fact that the amplitude satisfies an estimate like 
(HI]) ine. 

Now let us take a closer look at the critical points of ^ i— > '^{x,y,r],(^). We easily see that for 
X = y there is a unique critical for this function, namely = V- Combining this with (|45p and the 
implicit function theorem we see that for |x — y| sufficiently small there is still a unique complex 
critical point ^c, depending analytically on x, y, and rj. Further, T^gg(x, x, ^) = Id > and so 
it is still true that X^'^^(x, y, ry, ^c) > if |x — y[ < with sufficiently large (this is required 

to apply the complex method of stationary phase below). Making possibly larger still we can 
ensure that Aj^i and Cn are analytic and independent of on |x — y| < Thus we may apply 

the complex method of stationary phase (see [15], Theorem 2.8 and the remark after the theorem) 
to the rest of the integral in (|46p to obtain 



\x~y\<S/C^ 



ei^(-'y'V'^){AN{x, v) fj b^{x, e) + Cn(.x, v) </>) dx d^ 

(48) 



e^*(^'^'^'^=)(yl(x, (3; h) f' V{x, /3) + C{x, /3; h) cj)) dx + 0(e"^) 

x-j/|<5/Ci 

where we are using the notation /3 = {y,ri), and V{x,j3) is a new vector valued function. We will 
also write 

$(x,/3) = ^'(x,y,r/,^c). 

Expanding <I> about y = x we have 

^>(x, (3) = -ipxix, ■n){y -x) + {y - xf {-Tpxx{x, rj) + i'0^g(x, r])) {y - x) + 0{\y - x|^), 

and so if 5 is sufficiently small, by (|45p we have that X$(x,/3) > C^\y — xp for a positive constant 
on the domain of integration in (j48p when (3 is real. Now, combining (j46p . ()47p . and (j48p we 
have 

gi#(x,/3) (^^^^ ^. (a;, /3) + ^(x, /3; /i) 0) dx 

x-j/|<5/ci (^49) 

= O [Ne'l^ + {C^Nh)^^ , 

and by choosing N such that N < {eC'^h)~^ < + 1, which is possible when h is sufficiently small, 

_ 1 

we can change the right hand side of ()49p to 0{e 'c^). 
Now we change variables 

(x,y,r/) ^ {x,y,C) = {x,y,i^x{y,r])) 
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in (j49p (this is a vahd change of variables on a small enough neig hborhood of (0,0,^°) by (05])). 
Note that the new phase function, <I>(x,?/, ("), satisfies 

<^>^ix,x,C) =C, ^yix,x,C) = -C, $(x,x,C) = 0, 
and I^[x, y,Q > C \x — y\ for x near on the real domain. 

With this new phase, ()49p is now in an appropriate form to apply Sjostrand's definition of the 
analytic wavefront set (see [15j). Indeed, if there were were only a single amplitude and function in 
the integral (I49p we would be finished, as in However, since we have a covector field, and 
additional function, we must continue a little further. Here we view 0] as a system acted on 
by an operator with a vector- valued symbol {AV ,B). Note that we have 

ap(A?^,i?)(0,0,e°) = (5;v(0,e°)fr^(0,e°),Civ(0,e°)) = (u;(O,0o)eo,a(O,^o)) 

using the comments after (j43p and the fact that A(O,0o) = 1 since 6q is the velocity vector for a 
curve in T. Here Gp denotes the principal symbol. Now, following [18] we repeat the arguments 
above for every {xq,0) normal to {xq,SP) in a neighborhood W of (xo,^o)) which is possible since 
U is an open subset of Q, in order to obtain (0911 for different symbols and phase functions, with 
the new symbols satisfying 

ap{Aeh>e,Bem^,e) = (u;(0, 0) 0, a(0, 0)). (51) 

Now let {^1, ... , 9n-i} be a collection of vectors in W which form a basis for the space normal to 
. By the elliptic condition there must be some vector On = X]j=i l^'^^j ^ ^ such that 

n-l 

Glogu;(xo,e„) / J^^^Glogu;(xo,^j). (52) 

This provides us with n equations, 

f et*'("'^)(A,(x,/3;/i)/;6^(x,/3) + C7i(x,/5;/i)0)dx = ofe-^) , i = l, ...,n (53) 

J\x-y\<5/C^ ^ ' 

where for each z, <I>j is a phase function satisfying (|5Up . and (Aib^-jCi) satisfies (|5ip with 6 = 6i. 
We need one more equation to make an elliptic system, and for this we use the fact that is 
solenoidal on the interior of M. Thus 5/'^ = 0, and, as in |18j . if we take xo to be a smooth cuttoff 
function with support contained in M and equal to 1 near G M, then /iexp(^^>i(x, (3))xo^f^ = 0. 
Integrating this equality with respect to x yields 

= J /ie^*^(^''3)xo(x)(5r(x)dx = e^*^(^'^)(z($i)^(x,/3) +/iL»^(x))/;(x)dx. 

This is equivalent to 

J e:^*i(^'^)l)J(x, /3; cj)f^{x) dx = (54) 

where ap{D^){0,0,^^) = {(,^y ■ We now claim that the system of n + 1 equations obtained by 
combining (I53p and (I54p is elliptic at (0,0,^*^). By looking at principal symbols we see that this is 
equivalent to the condition that for a constant covector / = (/i, ... , /„) and constant <j) 

w{Q,9i){e,y fj + a{Q,9i)<l) = ^ = 1, ... ,n, and {fy fj = 

implies [/,</>] = 0. Indeed this is true by the ellipticity condition, ([52]) . and the fact that the 6i 
contain a basis for the space perpendicular to SP. Let us now rewrite the system in the more 
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compact form 

/ diag(e^*»(^''^))A(x,/3)[/^0](x)dx = O(e-^^) (55) 

J\x-y\<C 

where A(x, f3) is a matrix valued symbol acting on [f^, (j)] and C is a new constant. The ellipticity 
of the system at (0,0,^'') is simply equivalent to the invertibility of the principle symbol of A at 
this point. 

We would like to replace A in (j55p by the identity matrix so that, with this replacement, (j55p 
precisely shows, in the sense of [T5, Def. 6.1], that (xq, ^'^) ^ WFA{[f^, (p])- To do this we generalize 
the proof of proposition 6.2 from [15j to the case of operators with matrix valued symbols. Indeed, 
following [15] and [T8j we first define a system of ^'DO's in the complex domain by writing 

Op(A)[r,<A](y) = diag(e^(*'(^'^)-*ra)A(x,/3)[r,0](x)dxd/3. (56) 

These operators have differing phase functions, but by making appropriate changes of variables 
we can change them all to the same phase <1> and this does not change the principal symbol. 
Therefore, we may construct a parametrix for the system and following |15j use this parametrix to 
express Ide^^ as a superposition of functions Ae^^* modulo an exponentially decreasing function. 
Following now the same argument as is given for proposition 6.2 in [15j . but with matrix valued 
symbols, we have (possibly with a new constant C) 

/ ei'^^''^'^hd[f',(t,]{x)dx = Oie->^), (57) 

J\x-y\<C 

for P = {y,ri) in a neighborhood of (0,^*^). This proves that (0,^'^) is not in WFa([/*, i?^']) as 
originally claimed. □ 

Proof of TheoremUl Let / € L?'{M) be such that lG,x,wf = 0, and pick an intermediate manifold 
-^1/2 with M (£ Mi^2 <s Ml. Now extend / as zero outside of M and consider the decomposition 
pi|) on Ml/2 

We extend /l,/^^^ '^^1/2 ^^^^ outside of M1/2 and then apply Proposition H] to the pair 
[•/'mi/2' ^-^1/2] ^ L^(Mi) with M replaced by Mi/2- The proposition implies that the pair [flj^^^, 
is analytic in the interior of M1/2, and so by (I58p / is also analytic in the interior of Mi/2- Since 
/ is identically zero on M1/2 \ M this implies that / must be equal to zero on all of M. Therefore 
the kernel of Ig,x,w is just {0}, and so Ig,x,w is injective. 

□ 

4. Application to the attenuated ray transform of vector fields 

The Riemannian version of the attenuated ray transform of vector fields, see ([2]), ([3|), is given by 
the weighted transform (|4]) with a weight 

wAx,e) = g- -(^-''(^)'^-«W)<i^ (59) 

Compare this with ([3]), where a{x,6) is a given, usually non-negative, function on SM modeling 
the attenuation properties of the media, and the generator G is the generator of the geodesic flow. 
For the Riemannian case of the attenuated ray transform and its relation to the transport equation, 
we refer to [H]. In the analysis that follows, actually G does not need to be the geodesic generator. 
The weight (|59p satisfies the relation 

G log w = —a. 
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The results in section [2] then imply the following. 

Corollary 1. The elliptic condition is satisfied in any of the following cases: 
(a) w = Wa- is with a on Q. 

(h) U satisfies ^ and w = x^a, where Gx = 0, x\u 0; ^'^^ either a = a{x) is a function 
of X only, or a = aij{x)9''6^ ^ 0. 

In case (a), we integrate over all geodesies and we require the absorption a to be non-zero 
everywhere, although possibly anisotropic (dependent on 6 as well as x). In case (b), we restrict 
the ray transform to U using the smooth cut off x but we then require a to be isotropic and non- 
zero; or a non-degenerate quadratic form of 6 for any x. One can allow more general cr's in (b) but 
we only mention cases that might be useful in applications. 

Proof. If w is as in (a), then we can choose U = Q. Then Glogw = —a cannot be of the form ([9]) 
on for a fixed x because it never vanishes. 

Assume (b). Then either Glogw = —cr{x) or Glogw = —(Tij{x)9^9^ on the set X 7^ Oj ^-iid in 
particular in U . In either case, it is not of the form Q on fi [/. □ 

Acknowledgements. The authors thank J. Sjostrand for his help with certain aspects of the 
analytic microlocal theory presented in [15j . 
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